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The gravitational field of monopoles, cosmic strings and domain walls is studied in the quadratic 
gravitational theory R + aR 2 with a\R\ <C 1, and is compared with the result in Einstein's theory. 
The metric aquires modifications which correspond to a short range 'Newtonian' potential for gauge 
cosmic strings, gauge monopoles and domain walls and to a long range one for global monopoles and 
global cosmic strings. In this theory the corrections turn out to be attractive for all the defects. We 
explain, however, that the sign of these corrections in general depends on the particular higher order 
ffo^erivative theory and topological defect under consideration. The possible relevance of our results 
QtV the study of the evolution of topological defects in the early universe is pointed out. 

On 
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I. INTRODUCTION 

After the paradigm of the Hilbert's Lagrangian formulation of Einstein's theory of gravity it was clear 
hov^nne could consistently formulate other, higher derivative gravitational theories (that is theories in 
whi^^the field equations have higher than second metric derivatives). And such theories where indeed 
proposed and used as alternatives to Einstein's theory in attempts to unify other fields with gravity |l| 
and to remedy some of its seemingly undesirable consequences as, for example, at the classical level, 
thefuriavoidance of cosmological singularities JE| and, at the quantum level, the non renormalizability 
of t^^quantized version of general relativity 

OxH) of the main motivations for studying higher derivative theories comes from the semiclassical 
genetSl' relativity. There, it seems to be a matter of self-consistency to consider higher derivative terms 
in thje 1 gravitational Lagrangian since such terms arise generically in one-loop calculations 0] . Certainly 
this (Sif&ion of self-consistency is a delicate issue and, as Simon has recently suggested |y7 it needs to 
be recpnsidered if one wants to avoid undesirable semiclassical predictions such as unstable Minkowski 
spacetime. Another recent motivation for considering higher derivative gravitational theories is that 
suc^Jieories have arisen as low energy limit of several superstring theories Q . 

HitSier derivative theories are of interest to cosmology mainly because, even vacuum theories admit 
cosmological models which give rise to the, so called, Starobinsky inflation Q (see however Ref. || for 
a critisism on its consistency in the semiclassical limit), without fine tuning of the initial conditions 

In this paper we want to look at another topic of cosmological relevance namely, the effects of higher 
derivative theories on the gravitational field of topological defects as monopoles, cosmic strings and 
domain walls. These are objects that may have formed during phase transitions in the cooling down 
of the Early-Universe and may have played a key role in the formation of the large scale structure 
of the Universe mainly through their gravitational interactions Since their main interaction 

is gravitational, it is important to have an idea of what modifications one should expect in their 
gravitational field when the relevant gravitational theory has higher derivatives. Some work has recently 
been done in this direction, but only for gauge cosmic strings fTy] . This work shows that in the weak field 
limit only short range corrections to the Einstein theory arise which are associated with the presence 
of additional massive fields in the spectrum of higher derivative theories. However, this is not expected 
to be in general true, especially for global topological defects which are extended field configurations 
and not localized as the gauge cosmic strings. 

In this work we have in mind theories that can be separated in a part Co for the gravitational field g M „ 
constructed with geometrical scalars of the Ricci tensor R a b, and another part Cm containing matter 
fields with standard coupling to the gravitational field 
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C= ^-C G (Ra b )+C M (9ab). (1) 

Hereafter k := 8ttG where G is the gravitational constant. For theories of this type it has been noted 
that they can be recasted into an equivalent theory of Einstein gravity interacting with additional 
matter fields [(H^-p"!]. However, as it was stressed by Brans jl5| and we shall explain in the next section, 
this equivalence is in general only at a mathematical level and not at a physical one. Nevertheless, based 
on such an equivalent system, Whitt |0| was able to show that the black holes of general relativity are 
the only black hole solutions of R + /^theories (no hair theorem). 

For the discussion of this paper we will deal with theories that have as gravitational part the following, 
often appearing in the literature, Lagrangians 

£ G = V=g(R + aR 2 + pR»»RH (2) 

and 



C G = V^9F(R), (3) 

where a,/3 are some coupling constants, g := detf^, and R = g^R^v Finally the F in Eq. (||) is in 
principle an arbitrary function of the curvature scalar R. However, later on we will take F to differ only 
slightly from the Einstein value R, that is F = R + aR 2 with a|i?| <C 1. See Ref. [Q for a treatment of 
the F = R + aR 2 theory in vacuum and the Ref. together with references therein for generalizations 
to arbitrary F(R) in the presence of particular forms of matter. 

The structure of the paper is as follows. Section II contains a brief review of higher derivarive 
theories to the extent needed in this paper. The field equations for the theories in ([[]), (@) and (|l|),(|^), 
are written down, and their spectrum is explained. With the procedure that enables the recasting of 
these field equations into Einstein type ones we obtain the basic result that is used in the Sec. Ill 
for the comparison of the gravitational field of global monopoles, cosmic strings and domain walls in 
Enstein's theory, and in the quadratic R + aR 2 theory with a\R\ <C 1. Section III also contains at the 
beginning a brief introduction to the topological defects. Finally in Sec. IV we conclude with a brief 
summary and comments. 

Throughout this paper we use the conventions Ti — c — 1, metric signature ( — h H — h), Riemann tensor 
R a bcd ■= -d d T a bc + and Ricci tensor R ab := R c acb . 



II. THEORIES WITH HIGHER DERIVATIVES 



A. Field equations 



We shall give now the gravitational field equations for the higher derivative theories given by the Eqs. 
(|), (|) and Eqs. ([[]), (|). The field equations for g^u are obtained by varying the action corresponding 
to Eq. ([!]) with respect to g^ v and contain derivatives of the metric up to the fourth order. For the case 
of Cq of Eq. (g) they read 



(1 + 2aR)(R flu - -g^R) + |i? 2 <?^ 



{2a + ^)g^R, p 'P -{2a + (3)R,^ 



2' 

[3 



-2n SSm 



(4) 



Notice that the trace of this equation is an inhomogeneous massive Klein-Gordon equation for the 
curvature scalar R 

(6a + 2/3)R, p ' p - R = kT^ m K (5) 
Finally, the field equations for the theory (Q) and ([|) can be written as 
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F'G^ = «TW + ^(F - F'R - 2F!*) + F.^, (6) 

where F' — dF/dR and = i? M „ — \g^ v R is the Einstein tensor. The trace of this equation is 

3F' P + F'R-2F = kT w . (7) 



B. Spectrum of quadratic theories - 
Weak gravitational limit 



We would like to stress here the fact that quadratic theories do not contain only the usual massless 
(long-range) spin-2 graviton field but also, in general, two massive (short-range) fields with spin-0 and 
spin-2. 

This spectrum can be easily recognized in the case of Cg of Eq. (||) when one writes the field 
equations in the linearized weak field limit using a convenient gauge (coordinate system). Indeed 
following Teyssandier fl?|| we have that g^ v can be decomposed in the weak gravitational limit (where 
9nu = ri^ + h^ u with {h^l << 1) as 

9v.v = Vnv + hffi +XV fJ ,u + '>P l x V , (8) 



with the field equations 

UhfJ = -2/t(T^ - \T n ^ v ), 
(U-ml) X = --KT, m 2 := 6a + 2(3, 
{U-m\)4^ u = 2k(T^ u - -Trjftu), m^ 2 := -(3, 

and the gauge conditions 



(9) 



d a {hfj - V) = 0, 

(r b ~ ^V ab ),ab = 0. (10) 

Here indices are raised and lowered with the Minkowski metric tensor and the operator □ is the 
Minkowskian one. One recognizes in Eqs. (|§|)-([l0|) the usual Einstein contribution hffi, that is the 
graviton field which has 2 degrees of freedom. Then, a scalar field x with mass mo, which obviously has 
one degree of freedom and appears as an overall conformal factor (in the considered approximation). 
Finally the massive tensorial field V-V^ with mass mi which turns out to have five degrees of freedom 

(note that in contrast to h^J its components satisfy only one gauge condition) and thus it posseses the 
stucture of a massive spin-2 field. In order to keep the "mass" parameters mo, mi real we shall demand 
the no-tachyon constraint 

3a + (3>0, (3<0. (11) 

We leave for the next subsection the case of the Lagrangian in Eq. (|J) where we will go beyond the 
weak gravitational limit and we will see that the spectrum of this theory consists of a graviton and a 
massive interacting scalar field. 



C. Reformulation of quadratic theories 



Interestingly enough, besides g M „ there is an alternative canditate for the metric field of the spacetimc 
|T^JT^| , namely the 7 M „ which is the inverse of where 



dCg 

dR^ v ' 



(12) 
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and 7 := detj^. 

In particular for the Co of Eq. (0) we have 



-Tf" := (1 + 2ai?) 5 ^ + 2(3R a0 g a »g 



otfi Bv 



(13) 



Expressing the theory in terms of 7^ via a Legendre transformation, one can reduce the order of the 
derivatives that appear in the field equations from fourth to second. But what is also important is, 
that the resulting theory takes the form of Einstein gravity for the metric 7 Ml/ plus some additional 
massive fields interpreted as matter fields. Thus the equation can be considered as a non linear 
decomposition of g^ v in the physical spectum of the full theory: a spin-2 massless field j^ u , a scalar 
field that appears as a conformal factor and is a linear function of i?, and finally, a tensor field related 
to R ab with 5 de grees of freedom (R ab is symmetric satisfying the 4 contracted Bianchi equations and 
its trace is essentially the previously mentioned scalar field). 

We can be more explicit in the interesting case of the theory in Eq. ([!]) where Cq is given by Eq. (|^) . 
Here, only an additional scalar field appears since 7 M1/ and g^ are conformally related. Indeed Eq. (O) 
gives 
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— ■ '=gF'g lw 



which implies that 



Defining a scalar field tp (not to be confused with the tensorial field ip^ of Eqs. (||), (||) via 

F' = cxp (^[^^j , 
the field equations (^) are written in the system 7^ as 

where, hereafter, hats denote quantities with respect to the metric 7^ and 
The potential U(ip) is given by 



(14) 
(15) 

(16) 

(17) 
(18) 



U{^) = -F'-\RF'-F), 



(19) 



which can be written as a function of "0 alone in regions where Eq. (|16|) is invertible. Finally the scalar 
field ip satisfies the equation 



Dip = 



1/2 



F 



2F - RF' 



(20) 



which can be checked to be equivalent to the trace (^) of the initial field equations 
These field equations (17) and (pp|) follow from the Lagrangian 



c 



1 

2k 



(21) 



which shows that the quadratic theory ([!]) is, loosely speaking, equivalent to "Einstein's" gravitational 
theory for the metric 7 a b plus an interacting massive scalar field tp plus the "peculiar" (not anymore 
usual) matter fields of Cm (gab)- They are indeed peculiar if jab is considered as the metric of spacetime 
since then the dependence of Cm on g ao implies, via Eq. pa), a non standard interaction of the metric 
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field 7 a ;, and of the field tp with the matter fields of Cm- More on this issue of equivalence will be said 

at the end of this section. 

Case: F(R) g R + aR 2 with a \ R \ < 1. 

We will now consider the interesting case where F(R) can be expanded as a Taylor series around R = 
and deviates only slightly from Einstein's theory 

F(r) = R + aR 2 + 0(R 3 ), a\R\ < 1, a:=^-L_ n - (22) 



2 



Assuming that the 0(R 3 ) terms can be ignored in this expression then the field equations (|l7|) and ( pp| ) 
simplify considerably. Indeed, in this approximation Eqs. (|lq), (123) imply 



6 xl/2 



i> » 1-1 ai?, (23) 

while the leading term in the potential £/(■)/>) of Eq. ( [i"9| ) is, assuming the no-tachyon constraint a > 0, 
a mass term 

UW) = imgV 2 + 0(^ 3 ), = -L (24) 

2 6a 



The metric 7 MI/ and the field ?/> are obtained from the field equations ( |17| ) and (|20|) which, to lowest 
order in the approximation (B2J), read 

(Q-^«Q) V V^( 7 «6). (25) 

Making use of Eq. ( |l5| ) one can finally obtain the metric g^. Notice that, to the considered approxi- 
mation, it does not matter which metric we actually use in r'f. it is more convenient, however, from 
the technical point of view to use j^. 

Concluding we arrive at the following result: 
For a given matter source T^ff' (gab), the metric g^ v in the quadratic theory (]l|) and (|J) with F(R) ss 
R + aR 2 and a\R\ <C 1 is given by 

ff/w = [1 + xbW> (26) 

where j^v is the metric in the Einstein's theory with source T^ff \j a i,) while the field x satisfies the 
equation 

(□- m 2 )x = _« T (M) (7a6)) m 2 :=± (27) 

with the □ operator taken with respect to the jab metric. 

This result follows directly from Eqs. (|2^), (|l5|) using the variable \ related to the field ip via \ '■— 
— 2[k/6] 1 / 2, 0. Note that from ( ^3| ) follows aR = — x/2 and therefore the condition a\R\ <C 1 for our 
approximation is equivalent to |x/2| <C 1. Finally, let us notice that in the weak gravitational limit the 
Eqs. (^6|) and ( p7| ) are consistent with the j3 = limit of the linearized equations @ and (||). 



D. Some remarks 



Based on the decomposition (|T^) several authors 13 12JTJ,|T^,|T^] have dealt with the question of 
whether quadratic theories are equivalent to Einstein's theory plus some additional fields. It seems 
that this may well be true for vacuum theories. However, as was pointed out by Brans [pL5[ | , (see also 
Refs. [^0[^l)), a subtlety appears in the case where usual matter is present. The problem is that the 
equivalence principle, a basic guide that one may use in constructing theories coupled to gravity and 
in particular to Einstein's theory, cannot be valid in both the original and the reformulated theories. If 
it is valid in the original theory, then a test matter in Cm of the Lagrangian (0) will follow geodesies 
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of the spacetime with metric g^ but, in general, will fail to do the same in the spacetime with the 
metric 7 M „. In this sense we are not entitled to consider the reformulated theory as Einstein's theory in 
the presence of some interacting fields. In the case that one is philosophically inclined to consider the 
7^1/ as the physical metric, while the g^ as some sort of unifying field, then the equivalence principle 
should be implemented in the matter part Cm of Eq. ([!]) using the metric 7^ in the place of g^. Of 
course then, this Cm will be non standard with respect to g^. 

Whether or not nature chooses to couple usual matter universally only to a spin-2 field (as the 7^) 
and not to a more composite one (as the g^v) is far from being experimentally testable. Trying to find 
an answer one may, however, employ some criteria of principle, as positivity of energy pofl . In any case, 
the use of new variables, as those of Eq. ([l3]) and Eq. (|f5|), which turn out to simplify technically a 
physical problem, is undoubtfully very useful even if it is not clear whether one can attribute to these 
variables a foundamental character. 



III. TOPOLOGICAL DEFECTS 
IN HIGHER DERIVATIVE THEORIES 

Cosmological defects are formed during phase transitions in the evolving Early Universe whenever 
the symmetry group G of the relevant field theory breaks down to a subgroup H so that the vacuum 
manifold M — G/H has some non trivial homotopy group J22|. Such a symmetry breakdown at an 
energy scale rj can be realized, e.g., with an n-component scalar field <j)^ having a Mexican-hat type of 
potential 

X " 

n« = -T(E^ (0 -^) a - ( 28 ) 

i=l 

The homotopic structure of the vacuum manifold depends on the number n of components of the scalar 
field and, thus, we may have the formation of domain walls for n = 1, cosmic strings for n = 2, monopoles 
for n = 3. These defects are respectively surface-, line-, and point-like configurations. Sufficiently away 
from these configurations, at distances d> 6, the scalar field (f>^ approaches quickly its vacuum value 
J2i 4>^<t>( 1 ' ~ rj 2 . Here 8 is the width of the core of these defects, of the order of m^ 1 where = rj\f\ 
is the mass of the scalar field's massive mode. Typically, for symmetry breaking at grand unification 
scale, 6 ps lCr 30 cm and nrj 2 ps 10~ 6 . 

Depending on whether the symmetry that breaks down is a gauge (local) or a global one we have 
respectively the formation of gauge or global topological defects. In the case of gauge symmetry there 
exists a well defined core, with width 6, where most of the energy of the topological defect configuration 
is localized. On the other hand, for global topological defects the components of the respective stress- 
energy-momentum tensor have, outside the "core" , a relatively slow fall off due to the gradients of the 
Goldstone modes of the scalar field </>W. Thus, global defects are extended configurations. The reason 
for this difference between gauge and global defects is that in the case of gauge symmetry the presence 
of gauge fields can compensate the gradients of the scalar field. Finally, in the case of discrete symmetry 
breaking, which gives rise to domain walls, there are no Goldstone modes and thus domain walls are 
localized configurations. 

Based on the above properties, we will make in what follows the following approximations: 

(i) Gauge topological defects and domain walls will be considered in the zero core-thickness approx- 
imation and thus their stress-energy-momentum tensors will have components with appropriate 
Dirac 5-fuctions. 

(ii) For global defects, we will make the cr-model approximation where the scalar field is fixed to 
its asymptotic vacuum value everywhere outside the defect. This is a sensible approximation at 
distances from the defect sufficiently larger than the "core" width S. 

In the following subsections we will obtain the gravitational field of cosmic strings, monopoles and 
domain walls in the quadratic theory R + aR 2 with a\R\ <C 1. For this we will make use of the result 
of the previous section (see Eqs. (pq), (p^)), stating that the metric in the quadratic theory, cLS(Q), is 
conformally related with the metric in Einstein's theory, ds 2 E y 

ds 2 {Q) =(l + x)ds 2 E) , |x/2|«I, (29) 



G 



with x satisfying the massive Klein-Gordon equation ( J27| ) in the ds 2 ^ metric. A consequence of Eq. 
(29) is that there will be a modification of the "Newtonian" potential equal to %/2. We will have below 
the oportunity to study its nature and its range in the case of topological defects, be them localized or 
extended sources. 

In general we shall restrict our attention to sufficiently large distances, d, away from the core, (d> <5), 
but we will keep in mind that a proper treatment at short distances requires a proper model for the 
core of the defect itself. In this way we will be able to use the existing results in General Relativity for 
the gravitational field of cosmic strings, monopoles and domain walls which were obtained by making 
use of the above approximations in model Lagrangians with symmetry breaking potential of the form 

i§. 



A. Global monopoles 

The stress-energy-momentum tensor of a global monopole configuration, in regions far away from the 
core, can be approximated by p3| , p4j 

2 

T l= T r~~^i T e e = T£ = 0, (30) 



while the respective metric in Einstein's theory of gravity is (approximately) given by [|23|,|24| 

ds 2 {E) = -(1 - A)dt 2 + (1 - A)~ 1 dr 2 + dti 2 

dn 2 := r 2 (d0 2 + sin 2 Odtp 2 ), 
A := 8nGr] 2 = nrf. (31) 

This metric corresponds to a spacetime with a solid deficit angle: test particles are deflected by an angle 
7rA/2 irrespective of their velocity and their impact parameter. Here it should be added that a more 
careful treatment that takes into account the actual behaviour of the field at the monopole core, 
shows that the metric ([H]) gets modified by terms which at distances r 3> 6 = (^/\r|)~ 1 correspond 
effectively to a negative mass term M c g, that is e.g. gu ~ (1 — A — 2GAl c g/r). According to numerical 



analysis 24 M c g — 6n\/Xr]. Thus, besides the topological deflection caused by the solid deficit angle, 
test particles experience also a repulsive radial force —GM c s/r 2 away from the monopole. 

The metric in the quadratic theory is given by Eq. ( p9|) with \ satisfying Eq. (|27|). Looking for 
spherically symmetric solutions we find that this equation for x = x( r ) reads 



1 d 

r 2 dr 



2 d 
dr 



m >X{ r ) = 7 



3(1- A)r 2 ' 

rh 2 := ml/(l - A). (32) 



Making use of the Green function for this equation, 

G(r,r') = [e-""''sinh(mr)e(r' - r) 

mrr ^ 

+ e-" lr sinh(mr')e(r-r')], (33) 

where the step function O(z) := {0, 1, 1/2} for {z < 0, z > 0, z = 0} respectively, we can write down the 
solution for x(r) in terms of the Exponential-Integral (Ei) and Hyperbolic-Sine-Integral (shi) functions 
H as 

2A 1 _ _ — _ 

x( r ) = ~, rr^^PK - ™ r ) sinh(mr) — e -mr shi(mr)]. (34) 

3(1 — A) mr 

Checking numerically the behavior of this function we find that its contribution to the 'Newtonian' 
potential y/2 is an attractive one. In particular, using the asymptotic behavior of the Ei and shi 
functions |25[| we find that at large radial distances r — > oo 

X{r) M ~l (^' (35) 
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which implies a long range potential, exerting on test particles an attractive force — (2A/3m 2 l )r~ 3 . 
Comparing this force to the repulsive force due to the core of the monopole we see that the former falls 
off faster by one power of r and thus is negligible at very large distances. It overcomes, however, the 
effect of the latter at a distance r as m^~ 2 /(A<5) and, thus, it can be the dominant force within the region 
S <C r <C niQ 2 /(XS) which will exist provided that m^ 1 ^> S. 

Finally, let us note that the expression ( fii]) diverges as r — * 0. This is due to the form of the 
energy-momentum-tensor in Eq. ( |30| ) which is not valid at distances comparable to the core of the 
monopole. 



B. Gauge Monopoles 

A gauge monopole is a spherically symmetric configuration with mass M and a magnetic charge g. 
Its stress energy momentum tensor can be approximated by 

nt _ MS(r) ( 5 /4tt) 2 



TI = 



4-7T r 2 r 4 ' 

T r r =-T$ = -T* = - . (36) 

We consider the case where the metric outside the core of the monopole matches to a Reissner-Nordstrom 
one, (see Ref. p(l for a recent review and new results on the gravitational field of monopoles) 



ds 2 



, 2GM Gg z \ , 2 
1 + — ^ dt 2 

r 47rr z / 



2GM Gg 2 



r 



-l 



— + 4^) dr 2 + dn 2 . (37) 



Since the source for the x field is the trace of the stress-energy-momentum tensor, only the mass 
term in Eq. ( [36] ) will contribute. Furhermore, if we consider distances sufficiently far from the monopole 
r>i> GM , the equation for \ approximately reads 



LA 

r 2 dr 



2 ^ 
dr 



.. i . . kM 8{r) 
™l\x(r) = —^. 38 



Demanding finiteness at radial infinity, this equation has as solution the Yukawa fuction 

KM e- mor , s 

XW=-- . 39 

127T r 

Notice that the Newtonian potential of the monopole will be modified by the ammount x/2 correspond- 
ing to an attractive potential exponentially decreasing with an e-folding term characteristic of a massive 
scalar field with mass tuq. 

It worths remarking that the short range corrections of Eq. ( |39| ) apply also to the external metric 
of spherically symmetric mass distributions Jl7| such as neutron strars, giving thus rise to "fifth force" 
terms. However, when one deals with black holes, the no-hair theorem for R + R 2 theories [Of implies 



that corrections of the type (39) are absent 



C. Global Cosmic Strings 

As we explained in the introductory part of this section global cosmic strings are extended line 
configurations. The stress-energy-momentum tensor for a straight, static, cylindrically symmetric global 
string lying along the z-axis is approximately given for r 3> S by 

T^T^T^-T 6 ^-^. (40) 

The respective exact solution for the metric in Einstein's theory has been found in p7fl . However, it 
is quite complicated for the purpose of solving the equation d27n for the field \- Furthermore, besides 
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this technical problem, the spacetime of a global string has true spacetime singularities |2^,|2j|, a fact 
that demands carefull checking of the range of validity of the approximation (a\R\ <C 1) on which our 
treatment is based. Instead, we prefer to work here in the weak field limit where the equation for the 
X field is in Minkowski background metric. 

In the weak field limit of general relativity the metric of the global string reads |2q] 



l-4G7iln(- 



{-dt 2 + dz 2 ) + dr 2 



1 



8G/iln(- + c) 
o 



d9 2 



(41) 



Here [i := irq 2 , S is the core width and c is a constant of order unity that may partially take into account 
a global effect of the string core. Studying the motion of test particles it is seen that the static global 
string exerts a repulsive force 2G/z/r pSfl . It is interesting to explore how this force is modified in the 
quadratic theory that we are currently considering. 
The equation that x statisfies in the weak field limit is 



ld_ 

r dr 



dr 



Sirr 2 



(42) 



The Green function for the homogeneous part of this differential equation, with the boundary conditions 
of finitness at the origin and at infinity, is easily found to be 



G(r,r')= - K o {m r)I (m o r')e{r - r') 
- I (m r)K (m r')Q(r' - r) 



(43) 



where O is the step function. Thus the solution of Eq 

KfX 



can be written as 



X(r) = ~— K (m r) [ I a (m r ) — 
3tt L J s r' 

f°° dr' 
+ I (m r) / K (m r') — 



(44) 



where we have introduced a lower cutoff at r = 8 to cope effectively with the divergence that appears 
in the first integral if we let r — > 0. This divergence is only due to the approximate form of the 
stress-energy-momentum tensor which as we have already stressed is not valid near the core of the 
string. 

The leading term in an asymptotic expansion of Eq. (E3) at large radial distances is 



X(r) 



h'jt 



Z-nm^r 



2 r 2 ■ 



(45) 



From this expression we conclude that the additional 'Newtonian' potential in the quadratic theory 
implies at large distances an attractive force — (K^/3m§)r~ 3 . Due to the slower fall off of the original 
repulsive force, the total force on test particles remains repulsive at large distances from the string. At 
around r ~ 1 the total force is expected to change sign. 



D. Gauge Cosmic Strings 

Gauge cosmic strings are, in contrast to global ones, localized line configurations. The stress-energy- 
momentum tensor for a static, straight along the z-axis, gauge cosmic string with line energy density /i 
is 

^ = ^- 2,(1-^2) ^ T - T «=° ^ 
with corresponding metric in Einstein's theory Q 

ds 2 {E) = -dt 2 + dz 2 + dr 2 + (1 - Kfi/2) 2 r 2 de 2 . (47) 



9 



Here the polar coordinates r, 6 have the usual range. This spacetime is everywhere flat except along 
the z-axis where the string is located. As one goes around the string one notices an angle deficit. This 
topological property has the consequence that test particles which localy do not feel any gravitational 
forces are, however, deflected by the string. 

Let us now turn our attention to the field x- It satisfies the equation 



1 d 

r dr 



d_ 

dr 



-m^ X (r)= , K(i , W (48) 



which can be easily solved by demanding for the field x finiteness at infinity and correct discontinuity 
at the origin. The solution reads 

^ = - 3,(lZ,/2) K ^ < 49 > 

which can be easily checked that satisfies Eq. (Q) using the small argument asymptotic behavior of the 
modified Bessel function K (z) sw — ln(z/2). Finally notice that the field x decays exponentially fast 
since at large distances K (rriQr) « (7r/2m r) 1 / 2 exp(— m Q r). Very close to the string the expression 



(49) diverges logarithmically. Again, as the physical cosmic string has a finite core this divergence 
should not appear in a more proper treatment near the core. 

From Eq. ([l9|) follows that the 'Newtonian' potential of the cosmic string spacetime in an at- 
tractive, short range one. The respective force that the string will exert on test particles is 
— [K/umo/67r(l — K(J,/2)]Kx(mor). Because of the large distance exponential fall off behavior of 
Ki{rriQr) oc (mor) -1 / 2 exp(— m^r) it follows that this force is significant only close to the string up 



-l 



to distances r ~ m . 

In closing this subsection let us remark that the result obtained here is in agreement with the recent 
result of Linet and Teyssandier |ll[] in the weak field limit where Kjj, <C 1. These authors have also 
obtained the cosmic string metric in the weak field limit of the quadratic theory ([!]) and (||) which 
contains also the massive tensorial field tp^ of Eqs. (^)-((l^). In particular they find that the effect of 
this field on the Newtonian potential is a repulsive one with a range set by the inverse mass mi of the 
field Vv- 



E. Domain Walls 



The energy content and the gravitational field of domain wall configurations in Einstein's theory has 
been studied extensively in the literature, see Ref. []30| and references therein. We will consider here 
Vilenkin's vacuum plain domain wall solution discussed in pi] . For a domain wall with surface energy 
density <r, lying on the \z\ = plane, the stress-energy-momentum tensor is 

T* = T* = Tl = -aS(z), T* = 0, (50) 

while the respective domain wall spacetime is described by the metric pl| 

<k(E) = (! - v \Af [~dt 2 + e 2vt {dx 2 + dy 2 )] +dz 2 , 

v := 2irGa = kct/4. (51) 

Note that some of the metric components are time dependent (no static solutions can be found). Test 
particles in this spacetime are repelled with a proper acceleration v away from the wall pT[ | , a property 
that we may deduce just by looking at the 'Newtonian' potential term in the gu component of the 
metric ©. Finally we should mention that at \z\ = v 1 an event horizon appears (3lj. In what follows 
we restrict our attention to spacetime regions with \z\ < v^ 1 . 

Although the metric in (ph) is time dependent we can find, however, static solutions to the equation 
(|27[) for the field x depending only on \z\. For such solutions, equation (p7[) reduces to the ordinary 



differential equation 

mS) x (:) = W(i). (52) 
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The homogeneous part of this equation can be easily transformed into a Bessel differential equation for 
X where x(l z l) := x(^)/^ using the new variable z = v — \z\. In this way we find that the solution 
for x(z) is, in regions with z^O, a linear combination of the terms I\{jtiqz) j'z and Ki(moz) /z where 
K\, I\ denote modified Bessel functions. The coefficients of this solution are determined by demanding 
finitcness at the horizon z — 0, while on the domain wall, z = 0, the field \ should be continous and 
have the appropriate discontinuity in its first derivative which, according to Eq. (j52|), is [^x] 2 =o = 4za 
Thus we finally obtain 



m T , to , 



v v 



/^[l-^ID 



X (z) = -2 —h(-) \ 1 ' ■ (53) 



1-4*1 



It is easy to check that this implies an attractive and short range contribution to the Newtonian 
potential. This cannot overwhelm the original repulsive potential of the domain wall except very close 
to the wall for \z\ < mp 1 . This should be obvious in writing Eq. ( p3| ) in the sensible approximation 
Too jv 3> 1 and near the domain wall where \ takes its largest values 

X(z) ps -— exp(-TOo|^|). (54) 
to 

Away from the wall the field \x\ decreases exponentially and at the horizon x attains the small value 



IV. CONCLUSIONS AND REMARKS 

We have dealt in this paper mainly with the higher derivative theory ([!]), (||) with F = R+aR 2 in the 
approximation a\R\ <C 1. We showed that one can simplify the problem of solving the corresponding 
fourth order field equations using the conformal transformation ( [l5| ) which leaves us with the system 
of field equations (0) , (|2(i| ) having only second order derivatives. This is formally a system of Einstein 
type equations plus the field equations for a massive scalar field with mass to,q = 1/ (6a) interacting non 
minimally with gravity. We then found that in the approximation a\R\ <C 1 the gravitational fields in 
the R + aR 2 theory and in the Einstein theory are conformally related according to the Eqs. (p6|), fl27|). 
Using this result we looked in Sec. Ill for solutions representing the gravitational field of monopolcs, 
cosmic strings and domain walls 

For localized topological defect configurations as gauge strings, gauge monopoles and domain walls 
we have found short range, attractive corrections. Their range ~ 1/too is characteristic of the pressence 
of the massive field. 

For extended sources as global monopoles and global strings we have again found attractive corrections 
but with a long range. Their fall off rate depends on the corresponding stress-energy-momentum of these 
defect configurations. In particular, for distances r 3> to-q -1 is found that the attractive correction to the 
Newtonian potential is w kT/(6to§) where T is the trace of the corresponding stress-energy-momentum 
tensor. 

A more detailed investigation of the gravitational effects of topological defects in more general higher 
order derivative theories is in progress and we hope to present the results elsewhere. We can however 
already here reestablish the above given conclusions in the R + olR 2 theory and at the same time extend 
them to the case of the theory (|l|),((2|) which, for j3 ^ 0, contains also a massive tensorial field. This 
will be done by making use of the linearized field equations (||), (|9|). First observe that the Newtonian 
potential &n will have, besides the Einstein term also the contributions and <I>^f f "'' ) from the 

scalar field x an d massive field ip^ LV respectively 

<£ N = *(f> + + *JH = \[-h[f + x (55) 
For static spacetimes we have from Eqs. (0) and (|9|) that 

V 2 ^ ) = |(p + P 1 +P 2 +P 3 ), 

(V 2 - TO^ = %- Pl -P 2 -P 3 ), 
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(V 2 - m?)$»H = -|(2p + P 1 + P 2 + P 3 ) 

m 2 :=6a + 2/3, m^ 2 := -/?, (56) 

where p denotes the mass density and Pi, Pj, P3 the principal pressures of the matter. 

Sign of forces: Look at the r.h.s of these equations. For the topological defect configurations discussed 
in the present paper Pi is > for gauge monopoles, for gauge strings and global monopoles, and 

< for global strings and domain walls. On the other hand p—J2i P > 0> while — (2p + J2i Pi) < with 
equality holding for domain walls. Thus the Einstein contribution in Eq. ( |55| ) is attractive for gauge 
monopoles, zero for gauge strings and global monopoles, and repulsive for global strings and domain 
walls. The ^-contribution is always attractive, while the ^/^-contribution is in general repulsive except 
for domain walls where it is zero. 

Range of forces: From Eqs. (|5^) it is clear that the Einstein term provides in general a long range 
interaction due to an effective mass density p + £\ Pj. It gives, however, a zero effect for gauge cosmic 
strings and global monopoles. For the contributions of x and terms we have that: 

(i) If the stress energy momentum tensor vanishes (or falls off sufficiently rapidly) outside a localized 
source then at distances d from the source, these contributions are of short range cx exp(—md)/d p 
where m stands for the mass mo (or mi) and p is a parameter depending on the symmetry of the 
spacetime and is equal to 0, i and 1 for plain domain walls, strings and monopoles respectively. 

(ii) If the sources are not localized then for the x> contributions there are two characteristic 
regimes: 

(a) At distances d 3> max(m ~ 1 , mj -1 ) the mass terms dominate over the derivative terms in the two 
last equations of (|56|). Thus, asymptotically at large distances we have long range contributions 

$^ rj -n(p - SjPi)/(6mo) = KT/(6mg), as we found in this paper, and ps n(2p + 

£;P)/(3mf). 

(b) At distances d <C min(mQ~ , mj~ ) the derivative terms dominate over the mass terms. The 
interesting thing to note here is that at such distances the total 'Newtonian' potential in Eq. ( |55| ) 
satisfies V 2 ^^ ~ 0. This has implications for the differentiability of the spacetime metric and 
implies drastic changes in the singularity structure of gravity at short distances. For example, 
the gravitational potential of a point massive particle is finite at the origin in contrast to the 1/r 
Coulomb behavior in the Newtonian theory. 

These considerations are in agreement with the results of the previous sections and the results of 
p"l| ] for gauge cosmic strings. They may be particularly relevant to the study of the evolution of 
topological defects in the very early universe: (a) for structure formation scenarios based on global 
defects where the long range modifications of the quadratic theories may play an important role; (b) 
for collisions of cosmic strings where the drastic short range modifications may change significantly the 
predictions of these simulations for the evolution parameters of a string network. Thus it is interesting 
to study further topological defects and collisions of cosmic strings in quadratic gravitational theories 
and implement appropriate modifications in future numerical simulations. The outcome of such an 
investigation confronted with observation, may, among other things, allow one to put constraints on 
the mo, mi parameters of quadratic gravitational theories. 
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